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Stability  properties  of  hydromagnetic  waves  (shear  and  compressional  Alfven  waves)  in  spatially 
homogeneous  plasma  are  investigated  when  the  equilibrium  particle  velocity  distributions  in  both 
parallel  and  perpendicular  directions  (in  reference  to  the  ambient  magnetic  held)  are  modeled  by 
kappa  distributions.  Analysis  is  presented  for  the  limiting  cases  which 

solutions  of  the  dispersion  relations  are  analytically  tractable.  Here  ^f^(a=ej)  is  the  ratio  of  the 
wave  phase  speed  and  the  electron  (ion)  thermal  speed.  Both  low  and  high  (3  (=plasma  pressure/ 
magnetic  pressure)  plasmas  are  considered.  The  distinguishing  features  of  the  hydromagnetic  waves 
in  kappa  distribution  plasma  are  (1)  both  Landau  damping  and  transit-time  damping  rates  are  larger 
than  those  in  Maxwellian  plasma  because  of  the  enhaneed  high-cnergy  tail  of  the  kappa  distribution 
and  (2)  density  and  temperature  perturbations  in  response  to  the  electromagnetic  perturbations  are 
different  from  those  in  Maxwellian  plasma  when  <  1.  Moreover,  frequency  of  the  oscillatory 
stable  modes  (e.g.,  kinetic  shear  Alfven  wave)  and  excitation  condition  of  the  nonoscillatory  (zero 
frequency)  unstable  modes  (e.g.,  mirror  instability)  in  kappa  distribution  plasma  are  also  different 
from  those  in  Maxwellian  plasma.  Quantitative  estimates  of  the  differences  depend  on  the  specific 
choice  of  the  kappa  distribution.  For  simplicity  of  notations,  same  spectral  indices  /Cn  and  k  ^  have 
been  assumed  for  both  electron  and  ion  population.  However,  the  analysis  can  be  easily  generalized 
to  allow  for  different  values  of  the  spectral  indices  for  the  two  charged  populations. 

[DOl:  10.1063/1.3132629] 


I.  INTRODUCTION 

Low  frequency  (lower  than  the  ion  cyclotron  frequency) 
and  long  perpendicular  wavelength  (longer  than  the  ion  gy- 
roradius)  electromagnetic  waves,  often  referred  to  as  hydro- 
magnetic  waves,  are  often  observed  and/or  invoked  to  ex¬ 
plain  the  phenomena  in  both  space  and  laboratory  plasmas. 
These  waves  and  their  stability  properties  in  Maxwellian 
plasma  have  been  investigated  quite  extensively  for  many 
years  by  many  authors  and  a  good  discussion  on  them  can  be 
found  in  the  plasma  textbook  by  Stix.’  In  collisionless 
plasma,  however,  particle  velocity  distributions  can  often  de¬ 
part  from  being  Maxwellian.  For  example,  in  naturally  cx:- 
curring  plasma  such  as  plasma  in  the  planetary  magneto¬ 
spheres  and  in  the  solar  wind,  the  particle  velocity 
distributions  are  observed  to  have  non-Maxwellian  (power- 
law),  high-energy  tail.^  The  distribution  function  that  can 
better  model  such  particle  vekx:ity  distributions  is 
known  as  the  generalized  Lorentzian  or  the  kappa 
distribution^  with  functional  dependence  of  the  form  fo(v) 
1 For  finite  values  of  the  spectral  index 
K,  the  kappa  distribution  has  power-law  tail  at  velocities 
larger  than  the  thermal  velocity  0,  and  it  approaches  a  Max¬ 
wellian  distribution  [~exp(-i;^/ //^)]  in  the  limit  as 
Typical  values  of  k  for  space  plasmas  arc  in  the  range  of 
2-6.  In  the  last  several  years,  many  authors  have  studied 
electrostatic  and  electromagnetic  waves  in  spatially  homoge¬ 
neous  plasma  using  different  types  of  kappa  distributions  for 
the  equilibrium  state.^*'*'  In  a  recent  paper, wc  studied  low 
frequency  (lower  than  the  ion  cyclotron  frequency)  and  long 
perpendicular  wavelength  (longer  than  the  ion  gyroradius) 
electrostatic  waves  in  spatially  inhomogeneous,  current¬ 


carrying,  anisotropic  plasma,  where  the  equilibrium  particle 
velocity  distributions  were  modeled  by  different  kappa  dis¬ 
tributions.  In  the  present  paper,  we  investigate  the  stability 
properties  of  the  hydromagnetic  waves  in  spatially  homoge¬ 
neous  plasma,  where  the  equilibrium  particle  vekK'ity  distri¬ 
butions  in  both  parallel  and  perpendicular  directions  with 
respect  to  the  ambient  magnetic  field  are  mcxieled  by  kappa 
distributions. 

The  paper  is  organized  in  the  following  way.  In  Sec.  IL 
we  describe  the  general  mathematical  formalism  leading  to 
the  derivation  of  the  dispersion  relations  for  hydromagnetic 
waves  in  kappa  distribution  plasma.  In  Sec.  Ill,  we  analyze 
the  dispersion  relation  in  various  limits  and  describe  the  sta¬ 
bility  properties  of  the  well-known  hydromagnetic  waves. 
Both  low-/?  (/?<$!)  (y^?=plasma  pressure/magnelic  pressure) 
and  high-y3  (/3>  1)  plasmas  are  considered.  In  Sec.  IV,  we 
discuss  the  distinguishing  features  of  the  hydromagnetic 
waves  in  kappa  distribution  plasma  and  offer  some  physical 
interpretations. 

II.  MATHEMATICAL  FORMALISM 

We  adopt  a  Cartesian  coordinate  system  whose  c-axis  is 
along  the  ambient  uniform  magnetic  field  Bq  and  consider 
small  amplitude  electromagnetic  perturbations  represented 
by  (E ,B)cxp(-i(i)t ik  y-^ik^\z)  with  A:||>0  and  >0.  Elec¬ 
tromagnetic  modes  in  plasma  can  be  described  in  terms  of 
any  three  of  the  six  field  variables  (E,B)  by  eliminating  the 
other  three  with  the  help  of  k  X  E=(a)/c)B  and  k  B=0.  The 
three  field  variables  that  we  consider  to  be  physically  most 
meaningful  for  hydromagnetic  waves  are  and  B,.  The 
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magnitude  of  E.  is  a  measure  of  the  non-MHD  (magnetohy- 
drodynamie)  ebaraeicr  of  the  waves  and  corresponds  to 
the  compressional  component  of  the  magnetic  perturbation, 
which  characterizes  waves  in  finite-)^^  plasma.  Since  co<ck 
for  bydromagnetic  waves,  where  the  displace¬ 

ment  current  can  he  neglected  and,  consequently,  V  j  =  0 
(implying  quasineutrality)  can  be  assumed.  The  three  equa¬ 
tions  that  determine  and  are  then  the  quasineutral¬ 

ity  condition 

d\f„  =  o  (I) 

(t  a 


fMiv\sv;)  = 


r(/fii+ 1] 


xll+^ 


X 


+1) 


(4) 


The  function  normalized  to  the  equilibrium 

(unperturbed)  particle  density  Hq,  while  ^all  and  0^^  are  the 
parallel  and  perj^endicular  thermal  speeds,  respectively.  The 
parallel  and  perpendieular  pres.sures  in  the  equilibrium  state 
Pall  and  Pal  are 

r  2  /  2 /i’ll  \  1  ^ 

Pall  =  «o7’all='«„J  j-ff‘a"o(K,V 


and  the  perpendicular  (to  Bq)  components  of  Ampere’s  law 
given  by 


Here  is  the  eharge  and  is  the  perturbed  density  of  the 
charged  particle  species  a(  =£’,/),  and  j  are  the  components  of 
the  perturbed  current  density.  The  perturbed  quantities  arc 
calculated  from  the  perturbed  particle  distribution  function 
7«=7a(k,v,w),  as  shown  in  Eqs.  (l)-(3). 

We  refer  to  a  eylindrieal  coordinate  system  in  velocity 
space  with  its  z-axis  parallel  to  Bq,  so  that  V:^=v  ^  cos  cp,  Vy 
=  v  sin  (p,  and  1;^  =  v*Bo/Eo=V||,  where  (p  is  the  azimuth 
angle.  The  equilibrium  distribution  function  which  can 
he  constructed  from  the  constants  of  motion  (u^  and  Un)  of 
the  charged  particle  species  a,  is  taken  to  he  dcserihed  by  the 
produet  hi-Lorentzian-type  kappa  distribution  function  such 


Pa.  ->>0Ta^  =  ^'naj  (6) 

for  AC||>  1  /  2  and  k^>  ] .  Here  Ta^l  and  are  the  parallel 
and  perpendicular  effective  particle  temperatures,  respec¬ 
tively.  Two  other  parameters  related  to  and  which 
naturally  occur  in  the  study  of  electromagnetic  waves  arc 
Pa\\=^'rTPa\\^Bl=[2x^^/{2K^^-\)]{4^^mJ)QCFJBl)  and  /3„, 
^^TTpaJBl=[x  /{K,-\)]{47Tm^o(ii/Bl). 

In  our  choice  of  the  equilibrium  distribution  function  we 
have  allowed  the  possibility  of  different  values  of  the  spec¬ 
tral  index  in  parallel  and  perpendieular  direetions.  In  fact,  it 
may  he  reasonably  argued  that  /<■  >  /Cn  heeause  of  some 
equilibration  and  isotropization  in  the  perpendicular  plane 
while  preferential  acceleration  along  the  ambient  magnetic 
field.  Realistically,  the  spectral  indices  would  also  be  differ¬ 
ent  for  the  eleetron  and  the  ion  populations.  But,  for  simplic¬ 
ity  of  notations,  we  have  assumed  same  spectral  indices  for 
both  the  populations.  However,  the  analysis  pre.sented  in  this 
seetion  and  in  Sees.  Ill  and  IV  can  be  easily  generalized  to 
allow  for  different  values  of  the  spectral  indices  for  the  two 
populations. 

Solving  the  linearized  Vlasov  equation  by  standard  pro¬ 
cedure  (integration  along  the  unperturbed  particle  orbits  in 
tbc  ambient  magnetic  field  Bq),  we  first  obtain'^ 


7^(k,v,tu)  =  - 


Qa  ^  .^^iexp(-/X  eos  ip  +  inep)] 


^11  n--  +  2(co  -  /:||i;||)  ^  2 


w-  Ant'll  + 


<9  <9 

EzJniP-a){  ^11—  -  Y  I  + 


~  All  V  ,  ~ 

ck  k  c 


dv\  _ 


(7) 


where  7„(/z^)  is  the  Bessel  function  of  the  first  kind,  and  the  prime  notation  on  denotes  its  first  derivative 

with  respect  to  the  argument  /a^=A^i;  ^ /H^.  As  mentioned  earlier,  E^,  E^.,  and  By  have  been  eliminated  in  favor  of  E^,  and 
B^.  Using  the  identity 


+3C 

exp(-  ifji^  eos  <p)  =  2  (“  0%(m  Jexp(-  ipip)  (8) 


in  Eq.  (7)  and  earrying  out  the  cp- integrations  we  then  obtain 
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rlTT 

I  ‘ 


,  '  ,1  -2»n„-y)t 

n=-oc  -  Vll  +  I  \  ^^11 


nil,,  ku  u 

ck  k  c 


fl  d 

A'i|—  +  2(w-  /:||i;||)--Y 


(9) 


rlrr 

Jo 


d<f  cos  <pf„=—^  2 


(9  (9 


nil„  Lt  V 

ck  k  ,  r 


X 


(9  (9 

^■||~  +  2(w-A:||1;||)^ 
/9i;|,  r9i;  , 


f/oo(t^:-t^ii)- 


(10) 


.  -  27riq„^ 

J<psin<p/„= - —2. 

Jo  «=-= 


I  /l(l„  ~  J  r9 

^  - r; - — -2nn„—r\  + 

'nj-n  oj-kff„  +  nil„k  .  v  [  '  \  dvn  du ^ 


nil„  ^11  1 

- -jlB,+  i——J„J'B. 

ck '  k  c  "  "  ' 


A:||-  +  2(w-  A.||i;||)  2 

"^11  _ 


(II) 


which  arc  relevant  for  the  calculation  of  and  jy.  However,  for  the  description  of  the  low  frequency  and  long  parallel 

wavelength  mcxles  considered  here,  Eqs.  (9)-(ll)  can  he  simplified  by  taking  (cu~A||t;||+rtii^)"‘ 

=  (Aill^)“’[l-(a>-/:|,i;„)/;illJ  for  in  the  sum  and  using  the  identities  = 

n=0,  X„:^Qy„y^//i=0,  and  X„^q(7^)^/;i  =  0.  Keeping  only  the  leading  order  nonzero  terms  in  Eqs.  (9)-(l  1),  we  then  find 


27Tiq,A 


r^^fa= 

Jo 


E. 


2(1  -^5)-T-- 


Jo 

^11 


^  B,  All  V  ,  Vo/n 

(I  -Jl)—  +  i——  °  B. 


Sv\  6u-A|ii;||  ^V\\\  L  ^k  ^  k  c  (o-A||i;|| 


r?  r? 

Ai|—  -I-  2(tu-  Ai|i;||)— ^ 

ffV\\  Sv 


(/«o(l^!.t^ll)- 


(12) 


P 

Jo 


’  ,  7  277^/„I  /  £.  1  /  /' .u  n 

dip  cos  ipf„= - -)  I  - ; —  -  —  \J(ylo  +  i- - B, 

mji\\  I  \fu-A||i;||  ck  !  k^  c 


^Jo)  ^ii^ii 


-2{J\Y 
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A|1  +  2(tu-  AyDii)  2 

dv^\  Sv  ^ 


fAv^.  .i^ii). 


(1.^) 


r 


-  277/  11„ 

dip  sin  ipf^-—- - 

Am  a  V 


A 


^  ''0'  IV  ^  J  (y  Q 


k  ^  floJ 


d  d 

A’li—  +  2(w-  A:||t;||)— Y 

dvu  dv . 


ImP-v,)-  (14) 


The  expressions  for  n^,.  and  y\.  evaluated  from  Eqs.  (I2)-(14)  provide  adequate  description  of  hydromagnetic  waves  in  both 
{I5<  1)  and  high-/^  {p>  1)  plasmas  including  finite  gyroradius  effects. 

The  expressions  for  /T^,  and  jy  that  are  obtained  from  Eqs.  (I2)-(14)  (see  Appendix  A  for  some  details)  are 


n^,  2iq  f. 


A  ^4. 


!L 


(15) 


A  =  2  ^/a«07" 


A,{hA-A4(bA^Zl(U 


2<^^i 


C  —  to  B  y 
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«0  ^'ll  "0 
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ii»  <'„i 
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(16) 
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y’v  =  S  -  <S  ^-r  + 


a  ^fiK  ^'ll 


n 


201 


1  -A2{b„J  + 


2^11  <^ii 


A\{b„K) 


fin 


«  A.\,  «0 


(17) 


Here,  ^„=w/(/;||fi„||),  Z'J,^„)=dZJ,^„)ld^„,  b„^ 

=k\  /(20^),  and  the  Aj  are  defined  by 


<9 

|(finJ  =  2  -7o)^ 

Jo 


(\+V^/K^(t„J 


K  +1 


=  -l  1  -  ^b„^]b„^.  K^  >\, 


(18) 


i  +1 


ACi  3 


-  1 
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AC,  >  2, 


(19) 
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dv^Vjfy/f)-:^ 
Jo  '^I’l 
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2  AC  -  1 


(20) 
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(7o) 
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3  k' 

1  -  - - —b^Jb 


OK  I 


AC  >3 


(25) 


for  b„^<  1  (sec  Appendix  A  for  the  evaluation  of  the  inte¬ 
grals).  In  the  limit  as  ac^— 

which  is  a  notation  that  is  commonly  used  in  the  study  of 
Maxwellian  plasma,  and  Aj  go  over  to  the  corresponding 
exprcs.sions  for  Maxwellian  distribution.  For  example,  in  the 
limit  as  ac  — ^  _/^(/^^)^^xp(-/7j]  and 

^2ibaK)~^  ^o(bn)^^pi~bn)^  whcrc  Iq  is  the  modified  Bessel 
function.  The  function  Z^(^^)  is  the  plasma  dispersion  func¬ 
tion  when  the  parallel  velocity  distribution  is  given  hy  Eq. 
(4)  and  is  defined  hy 


16 
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(26) 
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- :(  1  - -b^^].  ac,  >2, 

AC  i  -  1  \  2  AC  -  2 


(21) 


and  by  the  analytic  continuation  of  Eq.  (26)  for  Im  ^^^0. 
This  is  similar  to  Z]^(^)  of  Summers  and  Thorne.^  The  series 
and  the  asymptotic  expansions  of  for  integer  values  of  aC|| 
are*^ 
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baK  I  ’  ^  1  ’ 


(22) 


Z,c(^J  =  I\'7r/(K||)(  1  - 
2kii  + 


«'ii+  1 


Nl 

U (27) 


6Kb„J  =  -^\  dv^v- - — - -j-— 


+  1 


AC  i  -  1 


1  - 


3ac^ 
— —t 

AC.  -2 


AC;  >2, 


Ai(b„J  dvAiUoY 

^a^Jo 


(23) 


dvi{\+vi/K,oi,r^-*' 


'C.  >2.  (24) 

AC  -  1  \  2  AC  -2 


z,iL)  = 


(l+^^'^ll) 

1, 


--1  1  + 


1 


L\  2k,- 


(28) 


where  /(ACy)  =  r(AC||-l- 1  )/[AC|]^^r(AC||+  1  /2)].  For  noninteger  (in¬ 
cluding  half-integer)  values  of  ACy,  the  asymptotic  expansion 
has  correction  terms, which  add  small  corrections  to  the 
real  frequency  of  the  excited  waves.  Wc  ignore  such  correc¬ 
tions  in  the  present  analysis.  The  linear  dispersion  relation  is 
obtained  hy  substituting  Eqs.  (15)-(17)  into  Eqs.  (l)-(3)  and 
demanding  nontrivial  solutions  of  E,,  and  B^. 


052106-5  Hydromagnetic  waves  and  instabilities  .. 


Phys.  Plasmas  16,  052106  (2009) 


III.  ANALYSIS  OF  HYDROMAGNETIC  WAVES 
A.  Shear  Alfven  waves 


In  the  study  of  shear  Alfven  waves,  eompressional  com¬ 
ponent  of  the  magnetic  perturbation  (5,)  may  be  ignored, 
while  the  magnetic  field  line  bending  represented  by  is 
more  important.  This  simplifies  the  analysis  considerably. 
We  assume  but  retain  finite  ion  gyroradius  (/?,><  1) 

effects.  After  substituting  the  appropriate  expressions  for  n„ 
and  into  Eqs.  (1)  and  (3),  we  obtain 


L 


()" 


2m,,  ^1, 


m,0l  ck 


(29) 


AAh.^)  co^ 
^11 


~  Ai{h,Jwck,~ 

n,  = 


(30) 


where 


+  — ^ — All 


AC  -  1  \  1  -  A2(/7j 


hi. 


A. 


(31 


Here  electron  terms  that  are  of  the  order  of  m^hrii  in  com- 
pari.son  with  the  ion  terms  have  been  neglected.  Also,  we 
have  used  the  definitions  of  b,^  and  V^  =  /?Q/(47rw,/io),  where 
is  the  Alfven  speed.  By  combining  Eq.  (29)  with  Eq.  (30) 
we  then  obtain  the  kinetic  dispersion  relation  for  shear  Al¬ 
fven  waves,  which  can  be  written  as 


1  -H 


^  \  ibjf(f 


^ A 


m,  b,^ 


+  g(b,,,fi)]^^k 


(32) 


We  now  solve  Eq.  (32)  for  different  limiting  values  of  |^„|. 
For  convenience,  we  shall  use  the  following  definitions  in 
our  presentation: 


EMa)  -  \  ^/( /<,!)[( +  I  «  I .  (33) 


LMa)  =  ^  7^/( /<„)[( /<ll  +  1  )//<l|]^.;(  1  +  I^J  I  • 

(34) 

First,  we  consider  |f,|>l  and  |^f,|>l,  i.e., 

<co/k^i  (cold  electrons  and  cold  ions).  Referring  to  Eq.  (28), 
we  realize  that  the  term  (Wt,/w,)(f^||/ ^?^,)A2(^/k)2^(^,)  in  Eq. 
(32)  may  be  neglected  in  comparison  with  as  mj rui 

1  and  typically.  Using  the  leading  terms  in  the 

asymptotic  expansion  of  Z^(^^),  we  then  find 


1  + 


k’lK  ^'ii  ^A 


+  ^(^>y3) 


-I- 


Ai{b,Jch\ 

u  2  .2i/2U 

biK  ^pe  ^11  ^A 


+  ^(/V.yS)]  =  0. 


(35) 


where  w^^^  =  4'rTe^nQl ni^.  Without  the  imaginary  term  that  ac¬ 
counts  for  wave-particle  interaction,  Eq.  (35)  yields 


klvl  A,(bJ\+{cVjwlX\+g(b,M 

Since  |fj 1  requires  [(2/C||- 1  )/2aC||]^^||<^//i,./aai, 
<^1,  [(2/C||“  1  )/2/C||]/3,||<^  1  and  since  we  can  assume 
Pnw  typical  plasma,  the  dispersion  relation  (36)  is  valid 
for  low-/3()8<^  1 )  plasma.  For  such  plasma  we  may  take  g 
1 .  Then,  Eq.  (36)  becomes 


_ i 

klvl~  A,{bJ\+(ch^/co;j' 


(37) 


When  AC  — Eq.  (37)  becomes  identical  to  the  dispersion 
relation  for  the  inertial  shear  Alfven  wave  that  was  derived 
earlier  for  Maxwellian  plasma. So,  Eq.  (37)  is  the  dis¬ 
persion  relation  for  inertial  shear  Alfven  wave  in  kappa  dis¬ 
tribution  plasma.  The  inertial  shear  Alfven  wave  has  been 
suggested  ^  as  a  mechanism  for  particle  acceleration  just 
above  the  auroral  ionosphere  where  the  plasma  is  cold 
enough  so  that  mf.  We  may  note  that  without  the 

(c^k^J  term,  the  dispersion  relation  follows  directly 
from  Eq.  (30)  if  we  .set  E,  equal  to  zero  (MHD  limit).  The 
(c^k^  /  term  arises  from  the  non-MHD  feature  (E,  ^0) 
and  denotes  the  kinetic  mexlifiealion  due  to  finite  electron 
skin  depth  (cico^^).  It  may  also  be  verified  that  b^^ 
<c"k\i  for  the  low-/^  plasma  defined  above.  Since 
A^/ -\  and  (1  -^2)//?,^=  ac^  /  (ac  -  I )  when  bj^<  1,  Eq. 
(37)  further  reduces  to 


kjvl  \+(cV\/wl)' 


(38) 


which  is  identical  to  the  dispersion  relation  for  Maxwellian 
plasma.  This  is  expected  because  under  the  conditions 
\^n\^  I,  b„^^0  and  with  the  omission  of  the  wave-particle 
interaction,  details  of  the  velocity  distribution  do  not  matter 
and  so  whether  the  distribution  is  kappa  or  Maxwellian  is 
unimportant.  With  the  inclusion  of  the  imaginary  term, 
which  is  small  compared  to  unity,  the  approximate  dispersion 
relation  for  low-)^  plasma  with  is 


k;vl 


c^k^  lit?  ~  , 

2/- - 


,  2;  2  /  2 
-Hc  k 


(39) 


It  exhibits  electron  Landau  damping  of  the  inertial  shear 
Alfven  wave  in  kappa  distribution  plasma. 

Second,  we  consider  and  i.c.,  ^,||<^^o/A.j 

^e\\  (hot  electrons  and  cold  ions).  For  simplicity  of  pre.sen- 
tation  we  ignore  the  terms  proportional  to  Z^(^,).  Tliis  elimi¬ 
nates  the  ion  sound  mcxle,  which  is  of  no  intere.st  here.  Using 
the  leading  terms  of  the  .series  expansion  for  Z^(^^J  in  Eq. 
(32)  we  then  find 
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,  i4|(/?,  ^)  0}  /  t  n\ 

'  2,,2 


h  k-^V^ 

^IK  '^11 

2ach  -  1  Ai(bi^}  ^  2  _2| 
2/Cii-i-l  bi, 


2/Cii 

1+/: — -LMe. 


2/<-||+  1 

k^.P%\+g{h^.,m  =  0-  (40) 


Here  p,=cj{\i  and  c^^=p,,||/(w,/io)  =  [/<-||/(2/<-||- l)](m,,/OT;)^ii 
Without  the  imaginary  term,  Eq.  (40)  yields 


kM  L  AM2k,+  \ 


bi^ 


[\+g(b^,,m-  (41) 


The  eondition  requires  [(2aC||- 1  )/2/C||])8^ii 

and,  as  before,  [(2/C||- 1  )/2/C||])0,||<^  1 .  Because  of 
these  conditions  we  may  neglect  )8,i|  and  (considering 
in  typical  plasma)  in  g  and  rewrite  Eq.  (41)  as 


b,^ 


2ku-  1 


2ac„+1 


I .  ^  .  ?srllk 


X 


/fL-lV4  2/f||+17',j^ 


1  +  -  (j8,.  -  ;3hi) 


(42) 


I  ^i(^k) 

"/K  ^li  1'A 


x(  1  + 


^2(^ik)  +  I 

^^11  ^1  ^|(^(k) 


0  ^  a  n  ,[1+«(^-k,;3)]  =  0. 

2/<-||+  1  ffr  A^ibi^) 

Without  the  imaginary  term,  Eq.  (45)  yields 


(45) 


1  - 


2/C||  ^|>\i(/7j 

2/C||  +  1  >42(^»k)  _ 


(46) 


The  dispersion  relation  (46)  is  valid  for  [(2/C||- 1  )/2/C||])t?^|| 
<tm^lmi  and  [(2/C||- 1  )/2/C||]y3,|il>  1 .  If  we  further  use  the 
small-/?, ^  expansion  of  Ai  and  A2,  then  Eq.  (46)  becomes 


1  -H 


2/C||-  1  r,n 

+  !/?,„ 


/c^  -  1  \4  2/C||  +  1  T, 


X 


(47) 


In  the  last  step  of  Eq.  (42)  we  have  used  the  small-/?, ex¬ 
pansion  of  /^i(/?,v)  [refer  to  Eq.  (18)]  and  the  relation  k^  p^ 
=  [Ki/(K^-\)]{T^/T,Jbi^. 

If  {mJm,)[2K^i/{2K„-])]</3,f<]  and 
then  Eq.  (42)  describes  the  kinetic  shear  Alfven  wave  in 
low-/0  kappa  distribution  plasma  consisting  of  hot  electrons 
and  cold  ions.  Unlike  the  ideal  MHD  shear  Alfven  wave,  the 
kinetic  shear  Alfven  wave  has  E,(?^0)  assex^iated  with  it  and 
incorporates  finite  ion  gyroradius  effect.  In  the  limit  as 
/C||,/c  — we  recover  the  dispersion  relation  for  the  kinetic 
Alfven  wave  in  Maxwellian  plasma,^^'^*^  which  has  been  in- 
voked  for  plasma  heating"  and  for  auroral  particle 
acceleration.^*^  With  the  imaginary  term  in  Eq.  (40)  included, 
the  approximate  dispersion  relation  is 


Comparison  with  Eq.  (42)  indicates  that  the  role  of  the  elec¬ 
tron  temperatures  is  now  played  by  the  ion  temperatures. 
When  Aii“  A  2,  Eq.  (47)  describes  another  type  of  kinetic 
shear  Alfven  wave  in  kappa  distribution  plasma  consisting  of 
cold  electrons  and  hot  ions.  On  the  other  hand,  if  Aii“A^ 
>2,  Eq.  (47)  describes  nonoseillatory,  purely  growing 
modes,  which  are  driven  unstable  by  the  ion  pressure  aniso¬ 
tropy.  This  is  similar  to  the  fire-hose  instability^*'^^^  condition 
[Eq.  (44)].  The  imaginary  term  in  Eq.  (45)  represents  the 
Landau  damping  (both  electron  and  ion)  of  the  modes. 

Finally,  we  consider  |^,j<^l  and  |^^|<^1,  i.e., 

>(x}/k^^  (hot  electrons  and  hot  ions).  Using  the  leading  terms 
of  the  scries  expansion  for  Z^(fJ  in  Eq.  (32),  we  obtain  the 
dispersion  relation 


^11 


biK  .  1,22 

-  -  +  p 

A^ibJ  2/C|,+  l 


2x11 

i - —LM. 

2/f„+l  ^  " 


(43) 


It  describes  electron  Landau  damping  of  the  kinetic  Alfven 
wave  in  low-/^J  kappa  distribution  plasma. 

For  high-y8^  >  1 )  plasma,  on  the  other  hand,  Eq.  (42) 
describes  nonoseillatory,  purely  growing  modes  when 

(44) 


The  instability  condition  is  the  same  as  that  for  the  fire-hose 
(garden-hose)  instability^^"^^  in  Maxwellian  plasma  with  an¬ 
isotropic  electron  pressures  (Aii^^t'i)*  T^^c  growth  rate  of 
the  instability  is  different  for  kappa  distribution  plasma  if  the 
finite  ion  gyroradius  effect  is  retained. 

Next,  we  consider  the  reverse  situation  |^,|<^  1,  ^ 

i.e.,  /?,ii  (cold  electrons  and  hot  ions).  Using  the 

series  expansion  of  Z^(^,)  and  the  asymptotic  expansion  of 
keeping  only  the  leading  terms  in  Eq.  (32),  we 

find 


^11 


.  .  '”.^11  .  .  .  2/C|| 

1  +  +  I 

2/f||  +  1 


m 


UQ  +  ^A^LM) 


'".^11 


b.^ 


2K„+\(b 


+  /- 


.  2x|| 


+  ^42L,(^,) 
m,(r. 


2/f||+  1 

When  /?,v‘^  E  Eq.  (48)  reduces  to 
1 


^11 


(l+g)  =  0.  (48) 


(49) 


where  and  The  dispersion  relation 

(49)  applies  to  high-y3  plasma  since  Va^ 

quires  [(2x||- 1  )/2x||]Aii^ ^nd  [(2x|i- 1  )/2x||]A|i^  1. 
It  indicates  instability  due  to  pressure  anisotropy  when 
A|-Ai^2.  The  instability  condition  is  similar  to  Eq.  (44) 
except  that  both  electron  and  ion  pressures  play  roles. 
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B.  Compressional  Alfven  waves 

For  the  compressional  AH'ven  waves,  we  may  neglect 
in  comparison  with  which  simplifies  the  analysis.  Wc  also 
consider  long  perpendicular  wavelength  modes  (^^,^<^1). 
Keeping  the  leading  terms  in  and  substitut¬ 

ing  the  resulting  expressions  Tor  and  into  Eqs.  (I)  and 
(2),  we  find 


fTlf,  &L\\ 
mitTu 


ieE. 


!L 

j/^o’ 


(50) 


(1) 


.2 


“*■  ,2  I 

A!  ^ 


-22^, 


z'M„) 


n-ll 


\B, 


k^. 


Ik^K  - 


k.B. 


11^0 


(51) 


where  electron  terms  that  are  of  the  order  of  in  com¬ 

parison  with  the  ion  terms  have  been  neglected  and 
1  has  been  assumed.  It  is  important  to  note  that  Eqs.  (50) 
and  (51)  are  valid  for  k’||>1/2  and  k  >2.  As  before,  we 
study  the  stability  properties  of  the  modes  for  dilTerent  lim¬ 
iting  values  of  |^^|. 

First,  we  consider  |^,|^1  and  i.e., 

<^a)/k^\  (cold  electrons  and  cold  ions).  Using  the  asymptotic 
expansions  of  Z^(^^),  keeping  only  the  leading  terms  by  not¬ 
ing  that  m^lrrii^  1  and  typically,  and  then  combin¬ 

ing  Eqs.  (50)  and  (51)  we  derive  the  dispersion  relation 


2 


2 


.2 


(o"  a:,,  a:  ,  .  m.  k^  kjvl 


-  r 


^  ■  r 


2/<-||-  \  K  ^-2  /3.,||  k' 


(52) 


A  possible  root  of  the  quartic  equation,  which  is  consistent 
with  |4|>l,  is  obtained  when  /3l^{mif  m^,)(k\f4k^) 
X(/:j[\/^/(t)^)  term  in  the  equation  can  be  neglected.  This,  in 
essence,  amounts  to  neglecting  the  term  associated  with  E. 
on  the  right-hand  side  of  Eq.  (51)  and  the  condition  for  it  is 
/3l ^  m,)[k‘^ /{klk\)].  When  the  condition  is  satisfied, 

the  dispersion  relation  is 

(i)^  k^  k 

1  -  -r-7  +  (/8  -  A)  2  +  ^ 

I.2y2  ^  1.2 


-  r 


f<\\ 


A.  A' 


2 


2/C|,-  1  /c^  -2  All  A 


-LMe)=0. 


(53) 


damping  (magnetic  analogue  of  Landau  damping)  in  kappa 
distribution  plasma. 

In  the  parameter  regimes  |^,|>1  and  1^1  i.e.,  f^ii 

(hot  electrons  and  cold  ions),  using  the  leading 
terms  of  the  appropriate  expansions  in  Eqs.  (50)  and  (51), 
and  neglecting  compared  to  terms  of  the  order  of 

unity  or  larger,  we  find 


2k«-  I  k}c'^ 


PiW 


ieE^ 


He.  All 


(54) 


i.z 

(O  A.|| 


2/<-||  +  I  /3; 


L( 

^  k^V^  -2  2/<-||- I  /3,,|| 


2/c,,  k^ 


K  -22k„  -  I  k^ 


-  All  All 


2  k,, 


k^\2K„+\  /  ,  A-iiA.  r 


X- 


leE. 


(55) 


By  neglecting  the  imaginary  terms  in  Eqs.  (54)  and  (55),  we 
obtain  the  dispersion  relation 


2k„-\  kA 
2  k,,  +  I  (x)^ 


n. 


k^ 

-rT+  -  A) 

k^v\  2k^  ^^ 


2/Cii+l  A 


■  22k,- I  a, 


Vll  '  J 


^  k\  2k, 1  /);, 


2k^2K,-  1  A,, 


(56) 


When  kllk^<2lli,  and  /)„>  I,  one  of  the  roots  of  Eq.  (56) 
corresponds  to  the  “slow”  magnetosonie  mode  with  or 
~[(2/<’|j- 1  )/(2/<’j|+ 1  )]A:j|r^,  and  the  root  is  approximately 
given  by 


klc] 


- 


2k„  + 


2Aii'<’.  -2  2/<-|| 


2k,-  I 


I+/3  , - 


pI'  X 


I  2k,,  + 


2/<-||+l\  "  /3^,|  /<•  -2  2/<-||- I 

The  mode  becomes  unstable  if  the  inequality 


(57) 


which  describes  the  damped  compressional  Alfven  mode  in 
kappa  distribution  plasma  including  finite  pressure  elTects.  In 
the  cold  plasma  limit  ^0),  it  reduces  to  (i)=kV^. 

The  real  part  of  Eq.  (53)  is  the  same  as  that  in  Maxwellian 
plasma.*  The  imaginary  term  represents  electron  transit-time 


^  ^  2m:||+1  a.  a.  ^  2(m:,  -2) 

K  -  1  2/<'||-  1  1  r  (3  ^  All  ^ 


(58) 


is  satisfied.  The  other  ro<it  corresponds  to  the  “fasf’  magne- 
tosonic  mode  o)'--kV^  and  is  approximately  given  by 


052106-8  B.  Basu 


Phys.  Plasmas  16,  052106  (2009) 


k^Vl  2Be,K,-22K,-\' 

The  instability  condition  for  this  mode  is 

2/C||+l  Bel  ^  2(/ci  -2) 

2/C||- 1  I  +  /c. 


(59) 


(60) 


In  the  limit  as  k:||,ac  — instability  conditions  (58)  and  (60) 
reduce  to  the  ones  that  were  found  for  Maxwellian  plasma, 
and  the  instabilities  were  referred  to  as  the  field-swelling 
instabilities.^^'^^  For  some  finite  values  of  /C||(>l/2)  and 
^■^(>2),  the  range  of  values  of  for  instability,  as  defined 
by  inequality  (58),  can  be  much  narrower  and  the  threshold 
value  of  for  instability,  as  suggested  by  inequality 

(60),  can  be  much  lower  than  those  for  Maxwellian  plasma. 
In  the  case  of  weak  pressure  anisotropy  so  that  the  plasma  is 
near  marginal  stability,  the  imaginary  terms  in  Eqs.  (54)  and 
(55),  which  take  into  account  Landau  damping  and  transit¬ 
time  damping  by  electrons  and  ions,  should  be  retained  in  the 
analysis. 

Next,  we  consider  the  reverse  situation,  i.e.,  |^,|<^1, 
|^^|l>l  (cold  electrons  and  hot  ions),  and  assume  that  0^ 
~0.  In  this  case,  Eq.  (50)  yields 


ieE.  AC,  2AC||-i-lm.  o/ 

_ - -  -L _ 11 _ c _ Z  ^  ^  (61) 

"*,■^•11^1  ”  1  2/c„  m, Bo  Bo' 

Therefore,  we  may  set  the  right-hand  side  of  Eq.  (51) 
equal  to  zero.  Moreover,  we  observe  that 
<{kllk^){dj,^/vlM(2K,-\)l2K^;\{kllk^)Bi,  as  Then 

the  dispersion  relation  that  follows  from  Eq.  (51)  is 

k^ 


Alfven  wave  due  to  ion  pressure  anisotropy  [see  Eq.  (47)  and 
the  discussion  following  it]  and  which  is  the  same  for  both 
kappa  and  Maxwellian  plasmas.  On  the  other  hand,  for 
nearly  perpendicular  propagation  (An//:^  1),  the  instability 

condition  is 


-2  2/f||- 

AC^  -  1  2aC||  -I-  1 


+  A-.), 


(65) 


which  is  the  condition  for  the  mirror  instability  in  kappa 
distribution  plasma.  In  the  limit  as  ac||,ac  — we  recover 
the  instability  condition  for  the  mirror  instability  in  Max¬ 
wellian  plasma. For  some  finite  values  of  aC||(>1/2)  and 
AC  ^(>2),  the  factor  on  the  right-hand  side  of  Eq.  (65)  can  be 
significantly  smaller  than  unity  and  thus  the  threshold  value 
of  /3-^/Aii  rnirror  instability  can  be  significantly  lower  in 
kappa  distribution  plasma. 

Einally,  we  consider  |^,|^1  and  l^t^l^l,  i.e., 

>  (jjlk\\  (hot  electrons  and  hot  ions).  If  the  imaginary  terms  in 
the  series  expansion  of  are  not  included  in  Eqs.  (50) 

and  (51),  we  have 


/  ieE,  _  K,  Idj,  OjU, 

\  m,djj  2(kj  -1)\ 


(66) 


o/  k; 

^,2^,2  +  2A,2(/^i-A) 


k^  „ 


I  - 


2k«  + 


I  0t 


k^-2  2/C||  (^11 


h 

Bo 


AC^  2/C||+  1 

Ac^  -  1  2aC||  k^  A:|j^Q  \  6^| 


(67) 


X 


^  ^^11  Ai  ^^11+  ^ 

-2  2ac,|-  1  !_  2aC|| 


(62) 


Here  we  have  used  (^  /  ^|=[(ac^  -  1)/ ac  ][2aC||/(2aC||- 1 )] 
^(Ai  /  All)-  Referring  to  the  definition  of  L^(^,)  we  solve  Eq. 
(62)  for  0}  and  find 


.  ^i|B,ii  k^-22k\^-\  Bn  k^ 

2v'77/(/C||)  -  1  '<11+  I  B]^ 


Jil  ,  + 


k^. 


+  ^2  \  '  +  - 


X.  -  I  2k\\  +  1 

K_-22K\^  -  I  ^,11 


(63) 


Instability  (Im  6j>0)  occurs  when 


Ax  ~  All  \  ^ 


|2  \  ^  A  L 


AC^  -  1  2aC||-H  1  /3-  \ 
AC,  -2  2aC||-  1  /3,||  / 


<0. 

(64) 


For  nearly  parallel  propagation  (/:„//:  ^^1),  the  instability 
condition  becomes  Aii”Ai  ^2,  which  is  same  as  the  condi¬ 
tion  for  the  fire-hose  (garden-hose)  instability  of  shear 


respectively.  For  the  convenient  case  of 

that  E^=0,  the  dispersion  relation  that  follows  from  Eq.  (67) 

is 


2 1/2 


1+ Jr(;8^-/3,i) 


2k^ 


+  ,  2  2  L  (  I  “ 


AC  -  1  2aC||  +  1 

AC^-2  2aC||-  1 


(68) 


when  ^  /(^ii  is  related  to  The  instability  condi¬ 

tion  can  be  written  as 


^  / “  Ai 
k^^ 


/f  , -2  2/C||-  I  ;8„|| 


<  0. 


(69) 


This  is  similar  to  Eq.  (64)  and  so  the  discussion  following 
Eq.  (64)  applies,  except  that  in  this  ease  both  the  hot  elec¬ 
trons  and  the  hot  ions  participate  in  the  instability  process.  In 
particular,  for  nearly  perpendicular  propagation  (^]|//:  <^1) 
and  for  An  Aih  ^he  instability  condition  is 


052106-9  Hydromagnetic  waves  and  inslabililies... 


Phys.  Plasmas  16,  052106  (2009) 


K  -2  2k\\  -  1 

,V 

Pill 


(70) 


where  (3^  =  The  instability  condition  (70)  is  evidently 

different  from  the  corresponding  condition  for  Maxwellian 
plasma,  which  is  obtained  from  Eq.  (70)  by  taking  ihe  limits 
When  the  imaginary  term  in  the  series  expansion 
of  Z^(^,y)  is  included,  the  dispersion  relation  is  modified  as 


K  K 


X 


K  ,  -\  2x11+  1 

-  2  2x11  -  1  /3^| 


2k\\ 

1  + 

2xii  +  1 


-ImZ  ’ 


(71) 

and  it  takes  into  account  the  electron  and  the  ion  transit-time 
damping. 

For  all  the  cases  considered  above,  the  Landau  damping 
and  the  transit-time  damping  rates  of  the  hydromagnetic 
waves  in  kappa  distribution  plasma  tan  be  significantly  dif¬ 
ferent  from  those  in  Maxwellian  plasma  depending  on  the 
allowed  values  of  /<^||(>1  /2)  and  k  (>2).  Further  discussion 
on  it  is  presented  in  Sec.  IV. 

IV.  DISCUSSION  AND  PHYSICAL  INTERPRETATION 

We  have  classified  the  hydromagnetic  waves  and  their 
stability  properties  in  kappa  distribution  plasma  according  to 
the  relative  magnitudes  of  the  wave  phase  speed  and  ihe 
electron  and  ion  thermal  speeds  (|^^|  and  |^,|)  and  considered 
only  the  limiting  cases  for  which  solutions  of  the  appropriate 
dispersion  relations  are  analytically  tractable.  Such  analytical 
approach  is  obviously  limited  in  scope  and  applicability;  but 
it  helps  to  understand  the  important  basic  features  of  the 
waves.  A  more  rigorous  analysis  requires  numerical  solution 
of  the  dispersion  relations.  Conclusions  reached  here  for  a 
specific  choice  of  the  kappa  distribution  function  would  re¬ 
main  qualitatively  unchanged  for  other  possible  choices  of 
the  kappa  distribution.^**^ 

The  well-known  results  for  the  hydromagnetic  waves  in 
Maxwellian  plasma  are  easily  recovered  from  the  results  re¬ 
ported  here  by  taking  the  limits  /C||,/c^  The  comparison 
shows  important  differences  between  the  kappa  distribution 
plasma  and  the  Maxwellian  plasma.  First,  the  magnitude  of 
the  resonant  wave-particle  interaction  in  the  suprathermal  re¬ 
gion  of  the  velcKity  space  (i;||=6a/A'||>  6^^||)  is  considerably 
larger  in  kappa  distribution  plasma  than  in  Maxwellian 
plasma.  This  is  evident  from  Fig.  1  where  the  behavior  of 
-Im  Z"  as  a  function  of  |ca/(A:||f?||)|  is  shown  for  two  types  of 
velocity  distribution.  The  reason  for  the  marked  difference  in 
the  magnitudes  of  -ImZ'  (which  is  proportional  to 
at  i?||  =  aj/A||)  when  |6a/{A:||^|)|  >  I  is  that  the  slope  of  the  kappa 
distribution  in  Un-space  decreases  according  to  a  power  law, 
whereas  the  slope  decrea.ses  exponentially  in  the  case  of  the 
Maxwellian  distribution.  Henee,  both  the  Landau  damping 
and  the  transit-time  damping  (magnetic  analogue  of  Landau 
damping)  of  the  waves  are  enhanced  and,  consequently,  the 
threshold  values  for  the  excitation  of  unstable  hydromagnetic 
waves  in  kappa  distribution  plasma  are  increased.  Other  im- 


FIG.  1.  Comparison  of  -Ini  Z'  vs  x.  where  .r=|a>/{A;  <? )].  Solid  curve  is  for 
Maxwellian  dislribution.  Dashed  curves  with  decreasing  lengths  of  segntenis 
are  for  kappa  distributions  with  /r  =2,  4,  and  6,  respectively. 


port  ant  differences  arc  the  following.  In  the  parameter  re¬ 
gimes  1^,1  >  1  and  1 ,  ( 1 )  the  frequency  of  the  oscillatory 
stable  mexie  (e.g.,  kinetic  shear  Alfven  wave)  in  kappa  dis¬ 
tribution  plasma  is  different  from  that  in  Maxwellian  plasma 
[see  Eq.  (42)]  and  (2)  the  instability  conditions  for  the 
nonoscillatory,  purely  growing  cotupressional  mcxles,  which 
are  excited  due  to  pressure  anisotropy,  arc  different  from 
those  in  Maxwellian  plasma  [see  Eqs.  (58)  and  (60)].  The 
condition  for  the  mirror  instability,  which  is  excited  in  the 
parameter  regimes  \^i\<  1  and  L  is  also  different  [see 

Eq.  (65)].  The  differences  can  be  attributed  to  the  fact  that 
the  density  and  temperature  perturbations  of  the  charged  par¬ 
ticle  species  a  in  kappa  distribution  plasma  are  different 
from  those  in  Maxwellian  plasma  when  L  This  inter¬ 

esting  aspect  is  discussed  below. 

When  and  neglecting  the  imaginary 

term  and  keeping  the  leading  term  in  the  rest  of  the  senes 
expansion  of  Z^(f^)  in  Eq.  (15),  we  find 

iq,,  2/Cj|  +  1  ~  /  2/Ci|  +  1  \ 

- ^ - 21-  ^  (72) 

"o  2x]|  -  1  ■'  \  2x||- 1  )  Bq 

after  relating  and  to  '^nW  and  7^  .  This  is  evidently 
different  from  the  corresponding  expression  for  Maxwellian 
plasma,  which  is  given  by  the  limit  of  Eq.  (72).  Simi¬ 

larly,  by  calculating  the  perturbed  parallel  pressure  /T^n  from 
Eq.  (12)  we  find 

PaW  _  n 
PaW 

which  is  clearly  the  same  for  both  kappa  and  Maxwellian 
plasmas.  That  Pa\\l P(y\\  should  be  the  same  for  any  choice  of 
equilibrium  distribution  is  easily  understoexi  if  we  reeogni/e 

that  it  also  follows  from  the  parallel  (to  Bq)  component  of  the 

28  ^9 

linearized  momentum  balance  equation 

0  =  -  (^|^„||  +  +  (A:||(/)„||  -/)„  (74) 

Equation  (74)  is  derived  from  the  linearized  Vlasov  equation 
by  taking  the  appropriate  velocity  moment  and  by  neglecting 
the  inertia  term  for  the  considered  low  frequency  waves.  We 


■■ 


(73) 
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then  notice  from  Eqs.  (72)  and  (73)  that  p a\\^ P ca\^  ^ J 
This  suggests  a  nonzero  temperature  perturbation  f^n  in 
kappa  distribution  plasma  given  hy 


Tc\\_  2 

^rtii  2  Kii  -  I  \  ^i|7’„ii  Bq 


(75) 


r+^  p 


r[(/+i)/2]r[,c„-(/-i)/2] 


r(K-,i+ 1) 


(A2) 


In  contrast,  p^ii/Paii=^a^^o  ^^id,  hence,  T^ii=0  in  Maxwellian 
plasma  under  the  same  low  frequency  conditions.  That 

7^11=0  in  Maxwellian  plasma  also  follows  from  Eq.  (75)  hy 
taking  the  limit  /cn— 

We  next  discuss  pai/pa^ ,  which  plays  role  in  determin¬ 
ing  the  excitation  conditions  of  the  field  swelling  and  mirror 
instabilities  [see  Eqs.  (60)  and  (65)].  Calculating  Pai^Pni 
from  Eq.  (12)  and  keeping  the  leading  order  terms  when 
1  and  we  find 

P»i  2/C||  +  1  -  _  /c  —  1  2/C||+  1 

Pal  ^l|7^a||2«’||-  I  ^  \  ~  2  2k-||-  I 

(76) 

Subtracting  Hq  wc  also  find 


for  k:||>  (/- 1  )/2,  where  in  the  last  step  wc  have  used  the 
standard  integral^’ 

r  r'  ria)r(b) 

The  second  integral  after  partial  integration  and  change  in 
integration  variable  is 

r-  dv,  ^8^  ^  r- 

j-cc  w-^||f^||  dvn  dwJ_o,  io)-k„Vi,) 

_] _ ^  ^ _ 

(A4) 


/Cj  2k„+1T„^  ^ 

K  -2  2k’|| -  I  7„||  _  B(j 


where  w/(^||^?„||).  It  is  related  to  Z’^{^„)=dZ^{^„)/d^„, 
(77)  where  Z^(^J  is  defined  hy  Eq.  (26).  The  needed  integrals  in 
space  are 


28  29 

The  corresponding  expressions  for  Maxwellian  plasma  * 
are  obtained  by  taking  the  /Cii, /c^ limits,  and  the  differ¬ 
ences  between  the  kappa  distribution  plasma  and  the  Max¬ 
wellian  plasma  are  evident.  Quantitative  estimates  of  the  dif¬ 
ferences  depend  on  the  specific  values  of  k’||(>1/2)  and 

/Cl(>2). 

In  Appendix  B,  we  present  a  more  convenient  form  of 
the  perturbed  distribution  function  that  may  he  used  in¬ 
stead  of  Eq.  (12)  to  obtain  and  p^i  given  hy  Eqs. 

(72),  (73),  and  (76). 
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APPENDIX  A:  DERIVATION  OF  EQS.  (15)-(17) 


In  obtaining  7\,  and  jy  from  Eqs.  (12)-(l4),  we  need 
to  evaluate  a  set  of  integrals  in  (i;||,i;x)  space.  When  Eq.  (4) 
is  used  for/^,  the  needed  integrals  in  i;||  space  are 


dvnv[gM  and 
J  — oc  J  —oc 


dUi\  dgM 
<-kM„  du«  ' 


(Al) 


where  1=0,2  and  ^oo('^ii)-[l ^t^er 
integrals  in  V\\  space  are  either  zero  due  to  symmetry  or  can 
be  reduced  to  the  first  integral  after  partial  integrations.  The 
first  integral  is  evaluated  with  the  change  in  integration  vari¬ 
able  as 


dv^ujlh^, 

Jo  Jo 


.  I.  jisdhoa 

dv^v^(\  -7o)— 


du^u\JoJih^,  du^u\JoJQ 
Jo  Jo 


dh 


oO 


(A5) 


f  dv  f  dvj^vliJof^^. 

Jo  Jo 

f  duj^u^JJ'iffiM,  [  du^uUj’ifj^, 

Jo  Jo  dv, 

where  h^{v^^)  =  [\ +v\/{k  01  When  the  scries 

expansions  of  Jq  and  7]  are  used  and  the  leading  few  terms 
arc  kept,  the  above  integrals  are  of  the  form  Jodv [\ 
-\-v\/(k  where  /=1,3,5,...  and  p=\,2.  It  is 

evaluated  by  changing  the  integration  variable  and  then  us¬ 
ing  the  standard  integral  [Eq.  (A3)]  and  is  given  by 


Jo 

=  {K^'^o„y*'r[{i+  i)/2] 

xr[AC,  +  p- (/+ l)/2]/[2r(Ac  +  p)]  (A6) 

for  /CjL  >  [(/+ 1)/2]-/?.  Evaluating  the  leading  terms  of  the 
integrals  in  Eq.  (A5)  using  Eq.  (A6),  we  obtain  the  Aj  given 
by  Eqs.  (18)-(25).  The  expressions  for  and  jy  given  by 
Eqs.  (15)-(17)  are  obtained  when  the  above  results  are  used. 
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APPENDIX  B:  A  MORE  CONVENIENT  FORM  OF 
WHEN  AND 

Here  we  construct  the  portion  of  the  perturbed  distribu¬ 
tion  function/^  that  is  most  relevant  to  the  calculation  of 
and  when  and  \o)/{ki\  ^al|)  I  ^  ^  •  (Landau 

damping  and  transit-time  damping  are  considered  unimpor¬ 
tant.)  Since  |cu/(lj  1 ,  we  argue  that  will  be  a  function  of 
tbe  magnetic  moment  pi=m„v\  1{2B),  which  is  a  valid  adia¬ 
batic  invariant,  and  the  total  particle  energy  in  tbe  wave  field 
e.  Keeping  in  mind  that  tbe  equilibrium  distribution  function 


is  a  kappa  distribution  of  the  form  given  by  Eq.  (4).  we 
c boose 


f„{fi.e)  =  C\  I  + 


-(x  +1) 


I  + 


2(e  -  fiBi)) 


(Bl) 


where  C  is  a  normalization  constant.  Using  fx 

"7  'y  ^  • 

where  eo=w^(i;^  -}-i;,|)/2  and  <t>  is  tbe  electrostatic  potential, 
and  expanding  in  a  Taylor  scries  around  /x=/Xo  and  k 
=  ^0  we  obtain 


J 


2(/C|,+  I) 


1  + 


(cq  “  /^o) 


-I 


Qa<f> 


/rtO 


-1 


/Cil  +  1  ^ 


1 


(^0  “  Mo^o) 


(B2) 


Here =  additional  terms  represent  tbe  perturbed  distribution  function/^.  Using  0  =  /E,/A:[|  for  the 

considered  low  frequency  waves  and  returning  to  the  velocity  space  variables,  wc  find 


2(/C||  +  1 ) 


1 


/C||  1  . 


K  -h  1 


1 


+  i;  Tf/c  )  «’||  6^1  1  +  _ 


B.  I 

(B3) 
"0  I 


f ao(^\  i^  same  as  that  given  by  Eq.  (4).  With  this  convenient  form  of/^,  the  necessary  velocity  integrations  can  be 

carried  out  more  easily  by  using  Eq.  (A3),  and  it  may  be  verified  that  /T^,  and  p^  calculated  from  Eq.  (B3)  are  identical 
to  those  given  in  Sec.  IV  from  a  more  complete  form  of 


'T.  H.  Slix,  The  Theorx  of  Plastna  Wciiv.v  (McGraw-Hill.  New  York,  1962). 
pp.  196  204. 

’d.  Summers  and  R.  M.  Thome,  Phys.  Fluids  B  3,  t835  (1991)  and  refer¬ 
ences  therein 

V  M.  Vasyliunas,  J  Geophys.  Res.,  [Space  Phys.]  73,  2839, 
DOl:  1 0. 1029/J  A073i009p()2839  ( 1 968). 

"^R.  M.  Thome  and  D.  Summers,  Phys.  Fluids  B  3,  2117  (1991). 

^D.  Summers  and  R.  M  Thome,  J.  Geophys.  Res.,  [Space  Phys  ]  97, 
16827,  1>OI:10.1029/92JA01664  (1992). 

^Z,  Meng.  R  M.  Thome,  and  D.  Summers.  J.  Plasma  Phys.  47,  445  (1992). 

^S.  Xue,  R.  M.  Thome,  and  D.  Summers.  J.  Geophys.  Res..  [Space  Phys.] 
98.  17475.  D01:10.1029/9.3JA00790  (1993). 

^D.  Summers,  S.  Xue,  and  R.  M.  Thome,  Phys.  Plasmas  1,  2012  (1994). 

Summers  and  R.  M.  Thorne.  J.  Plasma  Phys.  53,  293  (1995). 

•’’m.  a.  Hellberg  and  R.  L.  Mace,  Phys.  Plasmas  9,  1495  (2002). 

"R  L.  Mace  and  M.  A.  Hellberg.  Phys.  Plasmas  10.  21  (2003). 

‘‘R.  L.  Mace,  Phys.  Pla.smas  10,  2181  (2003). 

'^R.  L.  Mace.  Phys.  Plasmas  11,  507  (2004). 

‘"^R.  P.  Singhal  and  A.  K.  Tripathi,  Phys.  Plasmas  13,  012102  (2006). 

'^F.  Xiao,  Q.  Zhou,  H.  He,  H.  Zheng,  and  S.  Wang,  J.  Geophys.  Res.,  [Space 
Phys]  112.  A07219,  DOl  1 0.1 029/2006 J  AO  12050  (2007). 

"•B.  Ba.su,  Phys.  Plasmas  15,  042108  (2008). 

'^N  A.  Krall  and  A.  W  Trivelpiece,  Principles  of  Plasma  Physics 


(McGraw-Hill,  New  York.  1973).  pp.  402-405. 

'^R.  L.  Mace  and  M.  A  Hellberg.  Phys.  Plasmas  2.  2098  (1995). 

Summers.  R  M.  Thorne,  and  H.  Malsumolo.  Phys.  Plasmas  3.  2496 
(1996). 

~^'C.  K.  Goeriz  and  R.  W.  Boswell,  J.  Geophys.  Res.,  [Space  Phys.]  84. 

7239,  DOI:l0.l029/JA084iA12p07239  (1979). 

'‘C.  K.  Goenz,  Planet.  Space  Sci.  32.  1387  (1984). 

"R.  L.  Lysak  and  W.  Lotko,  J.  Geophys.  Res.,  [Space  Phys  ]  101,  5085. 

DOI;I0.1029/95JA03712  (1996). 

“^A.  Hasegawa  and  L.  Chen,  Phys.  Rev.  Lett.  35,  370  (1975). 

‘"^A.  Hasegawa,  J.  Geophys.  Res.,  [Space  Phys]  81,  5083,  DOl;  10. 1029/ 
JA081i028p05083  (1976). 

“^M.  N.  Rosenhluth,  Los  AlanK)s  Laboratory  Report  No.  LA-2030.  1956. 
‘^S.  Chandra.sekhar,  A.  N.  Kaufman,  and  K  M.  Watson,  Proc.  R.  Soc.  Lon 
don.  Ser.  A  245,  435  (1958). 

‘^A  A.  Vedenov  and  R.  Z.  Sagdeev,  Izv.  Akad.  Nauk.  SSR,  Ser.  Khi.  3.  278 
(1958). 

‘^B  Basu  and  B.  Coppi.  Phys.  Rev.  Lett.  48.  799  (1982). 

'‘'B  Basu  and  B.  Coppi.  Phys.  Fluids  27,  1187  (1984). 

'^’a.  Hasegawa.  Phys.  Fluids  12,  2642  (1969). 

'P.  J.  Davis,  in  Handbook  of  Mathematical  Tunettons  with  Formulas. 
Graphs,  and  Mathematical  Tables,  edited  by  M.  Abramowitz  and  1  A. 
Stegun  (Dover,  New  York,  1965),  p.  258. 


